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The standard model of electroweak interactions is minimally coupled to gravity and the 
response of the spherically symmetric solutions -the sphaleron and the bisphaleron- to gravity 
is emphasized. For a given value of the Higgs mass AIh, several branches of solutions exist 
which terminate into cusp-catastrophy at some (Afn-depending) critical value of the parameter 
a defined by the ratio of the vector-boson mass to the Planck mass. A given branch either 
bifurcates from another one at an intermediate value of a or persists in the limit a — > where 
it terminates into a flat sphaleron or bisphaleron or into a Bartnik-McKinnon solution. These 
bifurcation patterns are studied in some details. 

I. INTRODUCTION 

The coupling of non-abelian field theories to Einstein gravity constitutes natural extensions of the Einstein- 
Maxwell equations. One of the surprising issues of such possibilities has been the discovery of regular, finite 
energy classical solutions in the Einstein- Yang-Mills equations : the series of Bartnik-Mckinnon jl] solutions. 

On the other hand the Weinberg-Salam model theory, which couples the SU(2)(g)U(l) Yang-Mills theory to a 
doublet of complex Higgs fields emerges more and more as the theory of the unified weak and electromagnetic 
forces. This particular lagrangian is a member of the family of Yang-Mills-Higgs models among which the Georgi- 
Glashow model is another distinguished example. 

One interesting property of the Weinberg-Salam lagrangian is that it admits an unstable, finite energy classical 
solution called the sphaleron Q . Generally, sphaleron solutions are expected to play a role in the understanding 
of baryon non conserving phenomenon which are allowed to take place in the Weinberg-Salam model When 
the parameter determining the mass of the Higgs field increases, additional solutions, the bisphalerons, bifurcate 
from the sphaleron |^,||. This feature seems to be related to the underlying non- linear character of the classical 
equations and to the spontaneous breakdown of the symmetry ^-j^ . 

It is therefore natural to couple the Weinberg-Salam model to the Einstein-Hilbert gravitational lagrangian 
and to study the response to gravity of the classical solutions available in the fiat limit. This problem has been 
partly investigated in |^ and more recently in a report on the aspects of Einstein- Yang-Mills-Higgs equations 
|]lO| but it was not treated in details. In particular the response of the bisphaleron solution to gravity was 
not investigated. In this paper, we reconsider the classical equations of the Einstein- Weinberg-Salam model for 
spherically symmetric fields (we assume the Weinberg angle, related to the U(l) part of the gauge group, to be 
zero). Then the gravitating sphaleron and bisphaleron are constructed and assemble into branches of solutions 
which evolve with a, the ratio of the vector-boson mass to the Planck mass. The evolution of the different 
solutions as functions of a and Mh are studied in details. 



II. THE EQUATIONS 

We consider the gauge theory for an SU(2)-Higgs doublet minimally coupled to the Einstein-Hilbert gravita- 
tional lagrangian : 

L = V^[Lg + Lm] (1) 

with 

Lm = -\f;^F'^^'' + {D,^)^{D'^'i>) - ^(<i>t$ - (3) 
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and with the usual definitions for the fields strenghts and covariant derivatives 

F^. = d^^A^ - duA; + geabcA'^Al (4) 

D^$ = 9^$ + g(A><,)$ (5) 

The matter part Lm of this Lagrangian approximates the Weinberg-Salam model of electroweak interactions in 
the hmit of vanishing Weinberg angle 6w (i.e. the gauge group SU(2)xU(l) is restricted to SU(2)). 

Here we will study the classical, spherically symmetric solutions of the Lagrangian L. In this purpose, we 
employ the Schwarzschild-like coordinates for the metric 

ds^ = ^A^Ndt^ + N-^dr"^ + r'^(de'^ + sin^ Odcf^) , (6) 

and we introduce, as usual, the mass function m{r) by means of 

Then we use the standard spherically symmetric ansatz for the spatial components of the gauge field (fields are 
static and Aq = 0) and for the Higgs fields (the notations of are used) 

gr ' ' gr gr 

and 

$1 = , *2 = ^ [H{r) + iK{r){raa)\ . (9) 

It is well known that this ansatz for the matter fields is plagued with a residual gauge symmetry. Along with 
we fix this freedom by imposing the axial gauge 

xM, = ^ fc{r)^Q . (10) 

It is also convenient to introduce the dimensionless coordinate x and mass function /i defined by 

as well as the dimensionless coupling constants a, e 

where G is Newton's constant, v is the Higgs field expectation value, Mh is the Higgs boson mass and Mw is 
the gauge boson mass. If finite, the quantity /i(oo) defines the mass of the solution. 

With these ansatz and definitions, it can be checked after an algebra that the classical equations of the La- 
grangian (0) are equivalent to the equations derived from the following two-dimensional action 

J dt dx{-C) (13) 

C = A\^{N + xN' - I) + a^M] (14) 
where the prime means the derivative with respect to x and the quantity AI is defined by 
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M= NVi + V2 (15) 
Vi^ Ua? + {fsf + 2x\{H'f + {K'f) (16) 

^2- ^(/l + fl - 1)' + ^^\H^ + 1)2 (17) 

+ {H{fA - 1) + i^/s)' + (i^(/A + 1) - Hfef . (18) 
The classical equations then reduce to the following system of six non-linear differential equations : 

^' = a^M (19) 

A'^2Aa^-Vi (20) 

X 

{ANf^y = ^A(I^) (22) 
{x'ANH'y Ja(||) (23) 

{x'ANK'y = Ja(||) (24) 
It is important to remark that the equations are still invariant under the continuous global transformation 

fA + iJb (exp(2zl7))(/^ + z/b) 

H + iK — > {exp{in)){H + iK) (25) 

where is a real constant. Since the regularity of the functional V2 at x = clearly requires f\{0) + /^(O) — 1, 
we can fix the above symmetry by chosing /a(0) = 1,/b(0) = 0. This partly fix the boundary conditions which 
will be discussed more completely in the next section. Then, only a change of sign of the Higgs field can still be 
chosen arbitrarily. 

Let us close this section by discussing the limit a 0. From the definition (12) it is clear that the limit of 
vanishing a can be considered in different ways. If we keep v fixed and let G — *■ we obtain the flat limit (gravity 
decouples). The appropriate parameter which defines the classical energy of the solution is 

Ec = \imx^oc\fJ-{x) ■ (26) 

It restitutes the physically- meaningful classical energy of the flat (bi-)sphaleron In order to study the 

equations in the limit v and G fixed, then it is necessary to rescale the radial variable x and mass [i 
according to 

V=- , P=- 27 

a a 

and to set a — Q afterwards in the equations. We then obtain 



^ -2A-{{^Y + {^Y) (29) 



y dy dy 

^(Aivf )=a1/,(/1 + /|-i) , J = AB. (30) 

In particular, the degrees of freedom H, K related to the gauge field decouples and, setting /s = in the equations 
above, the Bartnik-McKinnon (BM) equations |jl[ are recovered. In passing note that we have not succeeded in 
constructing solutions of ( |30| ) which are not related to the BM solution by mean of (p5|). 
The corresponding energy is given by 

/i(x) 

Ebm = llTtly^ocPiy) = lini:r^oo (31) 
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III. BOUNDARY CONDITIONS 



The regularity of the solution at the origin, the finiteness of the mass /i(oo) and the requirement that the metric 
(^) approaches the Minkowski metric for x — s- cx) lead to definite boundary conditions (BC) for the six radial 
functions /i, A, /a, fs, H, K. As far as the metric functions are concerned we have to impose 

/i(0) = , A{oo) = 1 (32) 

For the matter field equations two sets of BC lead to regular and finite energy solutions. 0] 



1. The sphaleron BC 

The flat sphaleron (and also the gravitational one) has fsix) = H{x) = 0. The remaining functions have to 
obey 

/a(0) = 1 , K{0) = (33) 
/A(oo) = -f , if(oo) = l (34) 



2. The bisphaleron BC 

The flat bisphalerons are characterized by the four non-trivial radial functions and we will see in the next 
section that they are continuously deformed by gravity. Taking into account the fixing of the global symmetry 
(p5|), at the origin the functions have to obey 

/a(0) = 1 , /b(0) = (35) 
iJ'(O) = , K{Q) = (36) 

(37) 

In the limit x ^ oo, they have to approach constants in the following way 

\xm^^oc{fA{x) + ifB{x))= exp(2i7rg) (38) 
lim^^oo(-ff(a;) + iK{x)) = exp(i7rg) (39) 

So, the solutions of this type are characterized by a real constant q G [0, 1[; this parameter has to be determined 
numerically and depends on e and a. 



IV. THE SOLUTIONS 



We first describe the solutions for e = 0.5, a generic value of e for which the sphaleron is the unique solution 
of the flat equations. The flat sphaleron is there for a = 0, with an energy Ec ~ 3.64, and gets continuously 
deformed for a > 0. In particular the function N{x) develops a minimum which becomes deeper while a increases. 
This branch of gravitational sphalerons, let us call it Si{a)^ exists up to a critical value a = as ~ 0.3095 and no 
solution of this type exists for a > as. However there exists a second branch of solutions that we call S'u(a) for 
a £ [0, Qfs]. For fixed a the solution of the branch Su has a higher mass /is = m(oo) and a deeper minimum of 
N{x) than the corresponding one on the branch Si. This is illustrated on Fig.l; the indexes l,u in Si^u refer to 
the lower, upper branch when comparing the mass. 

In the limit a — > as, the solutions 5„ and Si converge to a common limit. It has 

Ai(as)w 0.290 , iV™„(as) « 0.513 (40) 



4 



The transition from the branch Si to the branch Su is completely smooth. 

In the limit a — > 0, the solution on Su tends to the first solution of the BM series. This is recovered by rescaling 
the radial variable x according to y = a;/Q; and taking the limit a = 0, as explained in the previous section. The 
BM solution is well known but, for completeness, we present its profile on Fig. 2. It has 

lim„^oo^«0.83 , iV„,„« 0.242 (41) 
a 

We have studied the gravitating sphaleron for a few values of e and determined the corresponding critical value 
as- These are presented in Fig. 3. The special value 0^(0) « 0.376 agrees with the result of Wc notice the 
rapid decrease of as for the low values of e. 

We next discuss the solutions for a value of e which allows both sphaleron and bisphaleron solutions to exist. 
This occurs for e > 72.0 and here we choose generically e = 100.0 for which the energies of the flat bisphaleron 
and sphaleron are given by 

Eh « 4.88 , Es « 4.93 (42) 

These two solutions get deformed by gravity when a > and develop two branches of solutions which we will 
denote Bi{a) (with mass ^b) and Si{a) (with mass fig). They exist respectively up to a = ~ 0.2247 and 
a = as ^ 0.2218. At these critical values, we have respectively 

Ai6(afc) « 0.2080 , iV„„« 0.464 (43) 
/i.(a,) « 0.2047 , iV„„« 0.463 (44) 

Again, these solutions are continued by upper branches which we denote respectively by -B„(a) and S'„(a) and 
which coincide with Bi{a) and with Si{a) respectively aX a — ai, and a = as- The two bisphaleron-solutions 
corresponding to Bi and i?„ for a = 0.2 are presented on Fig. 4. 

All along the two branches the energy of the gravitating bisphaleron stays slightly lower than the one of the 
sphaleron. Of course both quantities become equal when a = ac, this is illustrated by Fig. 5. This figure also 
clearly indicates that, at the critical value as (resp. ab), the energies of the two gravitating sphaleron Si,Su 
(resp. bisphaleron Bi,Bu) solutions form a cusp. When bisphalerons are present, four solutions are available on 
some interval of a and there are two cusps. 

Completely similarly to the case e = 0.5, the branch Su converges to the Bartnik-Mckinnon solution in the limit 
a — 0. The scenario with the branch Bu{a) is different: decreasing a from ab we observe that the different radial 
functions composing this solution uniformly approach their counterparts of sphaleron solution Su- This occurs 
for a ~ ac~ 0.185 as illustrated on Fig. 6. On this figure the quantities /b(oo) and H{Q) which characterize the 
bisphaleron are plotted for the different branches. The minimal value Nmin of the function N{x) is superposed 
on the figure. 

In view of these results we can say that the branch i?„ of solutions bifurcates from the branch 5^ at a w 0.185. 

In order to have a qualitative idea of how the gravitating (bi)-sphaleron solutions behave for higher values of 
e, we solved the equations for e = 800.0. There the gravitating sphaleron and bisphaleron exist respectively up 
to as « 0.209 and ab « 0.221. The branch Bu bifurcates from 5„ at a « 0.07 

This suggests that, when e increases, the branch Bu bifurcates from Su for lower values of a and that the 
gravitating bisphaleron exists on an interval which becomes relatively larger than the interval of existence of the 
sphaleron. 

In the limit of the non- linear sigma model e = cxi, it is known [§j4[|^ that the sphaleron is discontinuous at 
the origin unlike the bisphaleron which continues to exist as a regular solution. We expect that in this case the 
bisphaleron branch Bu will exist up to a = where it will approach the Bartnik-McKinnon solution. 

Let us finally say some words about the stability of the various solutions. For e < 72.0 the flat sphaleron 
possesses a single direction of instability (a negative mode) and for e > 72.0 the sphaleron has two directions of 
instability while the bisphaleron has one |^,|l3|. In this respect the sphaleron (resp. bisphaleron) is interpreted 
as the minimal energy barrier |^ bewteen topologically different vacua of the Weinberg-Salam model for e < 72.0 
(resp. e > 72.0). 

The shape of the plot of the energy, with the two branches terminating into a cusp is typical for catastrophe 
theory (see e.g. [Q) and the use of arguments based on Morse theory suggests many useful information about 
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the stability of the different solutions. For instance the number of negative modes for the solutions on the upper 
branch exceed by one unit the number of negative modes for the solutions on the lower branch. Such a reasonning 
was demonstrated to be correct e.g. in [ p^ . 

Using the same arguments in the present context indicates that the solutions on the branch Si of Fig. 1 have 
one direction of instability while the solutions of the branch Su (and therefore also the BM solution) have two. 

For the solutions of Fig. 5 we have one (resp. two) directions of instability for the solutions on Bi (resp. Bu) 
and two directions of instability for the solutions on Si . On the branch Su the number of negative modes is equal 
to two on the interval a € [0, ad (i.e. before the branch B^ has bifurcated) and equal to three for a G [Oc, as] 
(i.e. after the bifurcation). The number of instable modes at the approach of the BM solution is then equal to 
two, irrespectively of the parameter e. 

Obviously these deductions would need to be confirmed by more elaborated calculations. 



V. CONCLUSION 



The coupling of the Yang-Mills-Higgs equations to gravitation often leads to interesting new properties of the 
available gravitating classical solutions | p5| , p^ . 

The classical equations of the Weinberg-Salam model possess a rich pattern of bifurcations of bispahleron 
solutions |^,|| from the sphaleron solution Q . It is therefore natural to attempt to study the critical phenomenon 
which occur in the Einstein- Weinberg-Salam equations. In response to gravity, parametrized by the quantity a 
defined in (12), the equations exhibit another type of critical phenomenon : the occurence of two branches of 
solutions which terminate at a critical value of a; this was noted in pO[ | for Mh = but the phenomenon occurs 
for generic values of Mh ■ 

Here we were interested only in global solutions (i.e. the metric function N has no zero on [0, oo]), but another 
interesting feature of gravitationally deformed classical solutions (soliton or sphaleron) is the existence of black 
hole solutions ||l5| where the function N{r) develops an horizon at some finite value r = rh, i.e. N{rh) = 0. We 
guess that flat (bi-)sphalerons could also be deformed in this way and produce sphaleron black holes. 
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Figure captions 
Fig. 1. 

The mass of the gravitating sphaleron solutions (represented by fig and by iJ.s/ct) and the minimal value Nm 
of A'' are reported in function of a for e = 0.5. 
Fig. 2. 

The profiles of the functions N, A, /a of the solution Su as functions of y = a;/a for a = 0.01 and e = 0.5. 
Fig. 3. 

The critical value as a function of e for the low values of e. 
Fig. 4. 

The profiles of the functions N, A, f^, fs, H, K of the bisphalerons Bi (in dotted) and S„ (in solid) for a = 0.2 
and e = 100.0. 
Fig. 5. 

The masses fig , IJ-b of the gravitating sphaleron (solid) and bisphaleron (dotted) as functions of a in the region 
of the critical value. The branch iJ„ stops at a ~ 0.185, as indicated by the star. 
Fig. 6. 

The values of H{0), /^(cxj), Nm for the gravitating bisphaleron in function of a. The corresponding value of 
Nm for the sphaleron is represented by the dotted line. 
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Figure 2 
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Figure 4 
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